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Doping Dependence of Correlation Effects in Ki_^Fe2_j,Se2 Superconductor: 

LDA'+DMFT Investigation. 

LA. Nekrasoyi, N.S. Pavlov^, M.V. Sadovskiii-^ 
^Institute for Electrophysics, Russian Academy of Sciences, Ural Branch, 

Amundsen str. 106, Ekaterinburg, 620016, Russia 

Institute for Metal Physics, Russian Academy of Sciences, Ural Branch, 

S. Kovalevskaya str. 18, Ekaterinburg, 620990, Russia 

We present detailed LDA'+DMFT investigation of doping dependence of correlation effects in 
novel Ki-a;Fe2-i;Se2 superconductor. Calculations were performed at four different hole doping lev- 
els, starting from hypothetical stoichiometric composition with total number of electrons equal to 
29 per unit cell through 28 and 27.2 electrons towards the case of 26.52, which corresponds to chem- 
ical composition Ko.76Fei.72Se2 studied in recent ARPES experiments. In general case the increase 
of hole doping leads to quasiparticle bands in wide energy window ±2 eV around the Fermi level 
becoming more broadened by lifetime effects, while correlation induced compression of Fe-3d LDA' 
bandwidths stays almost the same and of the order of ~1.3 for all hole concentrations. However 
close to the Fermi level situation is more complicated. Here in the energy interval from -1.0 eV to 
0.4 eV the bare Fe-3d LDA' bands are compressed by significantly larger renormalization factors up 
to 5 with hole doping increase, while the value of Coulomb interaction remains the same. This fact 
manifests the increase of correlation effects with hole doping in Ki_a;Fe2-ySe2 system. Moreover 
in contrast to typical pnictides Ki_3;Fe2-ySe2 does not have well defined quasiparticle bands on 
the Fermi leves but pseudogap like "dark" region instead. We also find that with the growth of 
hole doping Fe-3d orbitals of various symmetries are affected by correlations in a different way in 
different parts of Brillouin zone. To illustrate this we determine quasiparticle mass renormalization 
factors and energy shifts, which transform the bare Fe-3d LDA' bands of various symmetries into 
LDA'-fDMFT quasiparticle bands. These renormalization factors effectively mimic more compli- 
cated energy dependent self-energy effects and can be used to analyze the available ARPES data. 

PACS numbers: 71.20.-b, 71.27.-|-a, 71.28.-|-d, 74.20. Fg, 74.25. Jb, 74.70.-b 



I. INTRODUCTION 



The discovery of iron based high-temperature super- 
conductors []| high-temperature superconductors stimu- 
lated quite intensive research work ^-ISl- Recently an- 
other class of high-Tc superconductors isostructural to 
122-family of iron pnictides was discovered - iron chalco- 
genides K^FeaSez 0, Cs^^FesSes ^ and (Tl,K)Fe^Se2 
[9| . Values of superconducting critical temperatures Tc 
are comparable for both families of pnictides and chalco- 
genides and are about 30-50K pi]| - [l3t . Further interest 
to these ehalcogenides was stimulated by experimental 
observation of rather nontrivial antiferromagnetic order- 
ing with very high Neel temperature of about 550K and 
Fe vacancies ordering in the same range of temperatures 
in K0.8Fe1.6Se2 (the so called 245 phase) [ij]. Despite 
intensive experimental work there is still no consensus 
on the composition of the phase, responsible for high-Tc 
superconductivity in these systems. The most common 
point of view is that KFe2Se2 (122 phase) is the par- 
ent compound for superconductivity, while 245 phase is 
insulating (Gl llSlllq . Some other phases in this system 
were also reported [13|- Below we shall concentrate on 
electronic structure calculations for the parent 122 phase 
with different levels of hole doping. 

Crystallographically pnictides AFe2As2, and ehalco- 
genides Fe(Se,Te), AFe2Se2 are quite similar with main 
structural motiff - layers of Fe(As,Se)4 tetrahedra. Re- 



cently reported theoretical electronic band structures of 
Fe(Se,Te) [H and AFezAsz [Hillj are found to be 
nearly identical to each other especially for those bands 
crossing Fermi level. AFe2As2 and AFe2Se2 compounds 
are simply isostructural. However LDA (local density ap- 
proximation) electronic structure of AFe2Se2 differs quite 
remarkably from that of AFe2As2 as was directly shown 
inRefs. :i5„22,23. 

From the very beginning of studies of iron based su- 
perconductors it was recognized, that the account of elec- 
tronic correlations on Fe sites is rather essential for cor- 
rect description of the physics of pnictide materials [2J- 
[131. To this end LDA+DMFT hybrid computational 
scheme [28| was employed. The main conclusion was 
that correlations lead to simple narrowing (compression) 
of LDA bandwidth by the factor of the order of 2 or 3. 
This observation agrees rather well with variety of angu- 
lar resolved photoemission spectroscopy (ARPES) data 
on AFe2As2 compounds [y. Fermi surface maps obtained 
from ARPES experiments for AFe2As2 are quite similar 
to those obtained from simple LDA calculations: there 
are 2 or 3 hole cylinders around F-pointin the Brillouin 
zone and 2 electron Fermi surface sheets around (tt, tt) 
point. 

Up to now there are only few LDA-I-DMFT papers de- 
voted to Fe ehalcogenides 29ll3G | . Recently we performed 
the investigation of electronic structure of hole doped 
iron chalcogenide Ko.76Fei.72Se2 in the normal phase 
[3l| . inspired by available ARPES data, for this system 



[3J|, especially those obtained in Ref. [35|, using both 
the standard LDA+DMFT and the novel LDA'+DMFT 
computational approach [33l33j . The results of our cal- 
culations agree rather well with the general picture of 
ARPES, obtained in Ref. [H, with LDA'+DMFT show- 
ing slightly better agreement. We showed that this iron 
chalcogenide is actually more strongly correlated in the 
sense of bandwidth renormalization (energy scale com- 
pression by the factor of about 5 in the energy interval 
±1.5 eV) around the Fermi level), than the typical iron 
pnictides (with compression factor of about 2 or 3 [y]), 
though the Coulomb interaction strength is almost the 
same in both families. Moreover the Ko.76Fei.72Se2 sys- 
tem demonstrates absence of well defined quasiparticle 
bands on the Fermi level in contract to pnictides. 

In this paper we continue our LDA'+DMFT study of 
system. We investigate the evolution of correlation ef- 
fects upon hole doping, performing LDA' calculations at 
four doping levels: hypothetical stoichiometric composi- 
tion with 29 valence electrons per unit cell through inter- 
mediate values of valence electrons 28 and 27.2 downto 
experimentally obtained composition Ko.reFei. 72862 with 
26.52 electrons per unit cell. We shall demonstrate that 
the actual doping dependence of correlation effects on 
the electronic structure in iron chalcogenides is appar- 
ently more complicated than in iron pnictides and does 
not reduce to the simple picture of universal bandwidth 
renormalization (compression). 

The paper is organized in a following manner. In 
Sec.|ll]we discuss crystallographic structure, methodolog- 
ical and computational details of LDA'+DMFT. Com- 
parative study of LDA' bands and LDA'+DMFT spectral 
function maps within wide and narrow energy intervals 
around the Fermi level, together with orbitally resolved 
densities of states are presented in Sec. IIIIl Finally we 
summarise our results in Sec. IIVI 



II. COMPUTATIONAL DETAILS 

The K2;Fe2Se2 system is isostructural to Bal22 pnic- 
tide (for the last one see Ref. [1^) with ideal body cen- 
tered tetragonal space group I4/mmm. The K2;Fe2Se2 
has a=3.9136A and c=14.0367A with K ions occupying 
2a, Fe - Ad and Se - 4e positions with 2;5e=0.3539 Q. 
This crystal structure was used in band structure calcu- 
lations for Ko.76Fei.72Se2 within the linearized muffin-tin 
orbitals method (LMTO) ^ using default settings ^. 

To take into account local Coulomb correlations we ap- 
ply here LDA'+DMFT [Hii approach which is a mod- 
ification of the weU known LDA+DMFT method (28} . 
LDA+DMFT Hamiltonian is usually written as 



H = H 



LDA 



^Hub _ jjDC 



(1) 



The general problem with LDA+DMFT is that some por- 
tion of local electron-electron interaction for presumed 
correlated d shells is already included in the standard 



LDA (iJ^°^). To avoid its double counting due to Hub- 
bard interaction //Hub ^^^^ ^^^ ^^ subtract from iJ^DA 
the so called double counting correction term H^^ . (Ex- 
plicit expressions for i/^DA ^^^ jijUuh ^^^ j^g found else- 
where ^.) The LDA'+DMFT approach is the new at- 
tempt to solve the double counting problem, which is due 
to the absence of universal expression for H^'^ , since 
there is no explicit microscopic (or diagrammatic) link 
between the model (Hubbard like) Hamiltonian approach 
and the standard LDA. 

To be short the main idea of the LDA'+DMFT is to 
exclude explicitly from the charge density the contribu- 
tion of the presumably strongly correlated d shells on 
the initial step of LDA calculations. Then this redefined 
charge density (for some chosen orbital basis (pi{r)) 



p'(r) ^ J2 l^'^WI 



(2) 



i¥'id 



is used to calculate the local exchange-correlation energy 
E^^^ in LDA and perform the self-consistent LDA band 
structure calculations for correlated bands at the LDA 
stage of LDA+DMFT. Further on the local d-d electron 
correlations are taken into account within the DMFT. All 
states not counted as strongly correlated are then treated 
with the fuU power of DFT/LDA and full p in Ei^^^. 

Once this LDA' calculations with redefined charge den- 
sity were done, what is left for the correlated states out of 
interaction on the LDA stage would be just the Hartree 
contribution which can be written in the fully localized 
limit form (FLL), and which is the most consistent defi- 
nition of double counting term here (other forms can be 
also used [33|): 

HpLL = i^Und{nd - 1) - -^J^ndairida - 1), (3) 

where nda = J2,n "^Udyna = J2mi'^iidrn<j) IS the total num- 
ber of electrons on strongly interacting orbitals and num- 
ber of electrons per spin, Ud = J2a ''^d.a-- 

Effective five orbital impurity problem for 
Ki_i;Fe2-ySe2 within DMFT was solved by Hirsh-Fye 
Quantum Monte-Carlo algorithm j37j . at temperature 
280K. LDA'+DMFT densities of states and spectral 
functions were obtained as discussed in Ref. j33 |. 
Coulomb parameters were taken to be U—3.75 eV and 
J— 0.6 eV respectively [33], which are very close to 
calculated ones [33|. To define DMFT lattice problem 
we used the full (i.e. without any downfolding or 
projecting) LDA Hamiltonian, which included all Fe-3d, 
Se-4p and K-4s states. 



III. RESULTS AND DISCUSSION 

In Fig. 1 there we present the comparison of 
LDA'+DMFT calculated spectral function maps in the 
wide energy window ±2 eV along high symmetry direc- 
tions in the first Brillouin zone with the renormalized 



LDA' bands (grey lines) for Ki_j;Fc2_ySc2 at different 
hole doping levels n^. Renormalization (bandwidth com- 
pression) factor of LDA' bands here is only 1.3, so that 
bandwidth renormalization due to correlations is rather 
weak. The lower panel of Fig. 1 LDA'+DMFT shows 
data for stoichiometric KFe2Se2 compound with total 
number of valence electrons ne=29. For this composi- 
tion all quasiparticlc bands are rather well defined for 
this wide energy range around the Fermi level. For hole 
doped cases with ne=28, 27.2 and 26.52 (2nd, 3rd and 
4th panels from bottom correspondingly) we see that Fe- 
3d bands obtained from LDA'-I-DMFT become less pro- 
nounced with hole doping. Overall rigid shift of LDA' 
bands from stoichiometric case down to the most hole 
doped one is about 0.3 eV. 

In Fig. 2 we show LDA'-fDMFT spectral function 
maps in the vicinity (-0.5 - 0.2 eV) of the Fermi level 
with dominant orbital character of quasiparticlc bands 
denoted explicitly by black squares for xz,yz, black cir- 
cles for xy, white circles for 3z^ — r^ and white squares 
for x^ ~ y^ . We can see thet the orbital characters and 
forms of quasiparticlc bands in this energy interval (-0.1 
- 0.1 eV) change with the increase of hole doping, al- 
though characters of quasiparticlc bands located outside 
this region remain the same. Main orbital character of 
bands crossing the Fermi level is xz, yz and xy. Also from 
Fig. 2 one can conclude that close enough to the Fermi 
level (for all hole dopings) there are uniformly no well 
defined quasiparticlc bands (though some low intensity 
maxima of spectral density can still be seen). This fact 
tells that the Ki_a;Fe2-ySe2 for all hole dopings is more 
correlated than 122 pnictide system. 

This "pseudogap" like behaviour can be explicitly ob- 
served in Fig. 3 for all Fe-3d orbitals, where orbitally re- 
solved bare LDA' and LDA'-fDMFT densities of states 
(DOS) for all four hole doping levels are presented. Also 
by watching these DOSes one can actually see, that upon 
hole doping correlation effects become stronger. This fact 
manifests itself in a different way for orbitals of various 
symmetry. First of all, for all Fe-3d orbitals we observe 
narrowing of DOSes. For 3z^ — r^ (third panel from top) 
and x^ — y^ (upper panel) this narrowing is most evident 
at ±1 eV and ±0.5 correspondingly. For xz,yz (second 
panel from top) and xy (bottom panel) increase of nar- 
rowing with doping is mostly concentrated within the 
energy interval ±0.4 eV. 

To get more insight into the LDA'+DMFT self-energy 
effects on bare LDA' bands we have determined energy 
scale renormalizations and energy shifts for variety of sep- 
arate dispersions of bare LDA' band structure depicted 
in Fig. 4, which rather accurately fit bare bands to those 
plotted in the Fig. 2. Also in Fig. 4 we show the stan- 
dard LDA bands (dashed lines), just to emphasize that 
LDA' dis per sions are quite close to LDA ones (see also 
Refs. [3l| - (33| ). The obtained energy scale renormaliza- 
tion (bandwidth compression) factors and energy shifts 
(shown in brackets) are collected in the Table 1 for all 
hole doping levels. These results show more complicated 



picture of bare LDA' dispersion transformations, than 
those obtained in Ref. [31[ , where we have proposed that 
all Fe-3d band dispersions for the Ko.76Fei.72Se2 should 
be just compressed by a factor of 5 to get reasonable 
agreement with experiment. In contrast to our previ- 
ous work Ref. 3l[ here we fit LDA' bare bands exactly 
to maxima positions of the spectral function (see stars 
on Fig. 2). Thus wc get more detailed picture which in 
general agrees with our early conclusions of Ref. [31| . 

Here wc actually see, that different parts of bands in 
the first Brillouin zone acquire different (bandwidth or 
mass) renormalization factors, which change with dop- 
ing. For example the xy,xz,yz band (see Fig. 4) almost 
docs not have doping dependence at all (see Table 1). 
Renormalization for the part of xz,yz band (2) only 
slightly depend on doping. However, the xz,yz band (1) 
becomes monotonically more correlated (renormalization 
factor grows up to 4.2) while hole doping increase. Cor- 
relation renormalization of the xy orbital demonstrates 
nonmonotonic behaviour and is more pronounced for in- 
termediate dopings. The renormalization factor of the 
3z^ — r^ band near the F point abruptly increase up to 
4.7 at ne=26.52. Here one should note that all four dop- 
ing levels were treated with the same Coulomb interac- 
tion parameters. Despite this fact upon doping several 
orbitals of Ki_j;Fc2-i,Sc2 become more correlated (nar- 
rowed). This phenomenon is apparently related to the 
change of the correlated orbitals occupancy. The data 
collected in the Table 1 might be helpful for interpreta- 
tion of APRPES spectra in simple terms of renormalized 
bare LDA band dispersions. 



IV. CONCLUSIONS 

In this paper wc have performed detailed 
LDA'+DMFT study of correlation effects in 
Ki_2,Fe2-j/Se2 system at four hole doping levels: 
from hypothetical stoichiometric composition with 29 
valence electrons per unit cell through intermediate 
values of valence electrons 28 and 27.2 towards compo- 
sition Ko.76Fci.72Se2 with 26.52 electrons per unit cell, 
for which there are available ARPES data on electronic 
dispersions |35| . 

Within DMFT correlation effects are concentrated in 
the self-energy which provides two types of modifications 
of the bare spectra at each energy: broadening (lifetime 
effects) by imaginary part of self-energy and energy shift 
due to its real part. We have shown, that in rather 
wide energy window ±2 cV (in terms of LDA energy 
scale) for all dopings the lifetime effects are relatively 
weak and renormalization (compression) of quasiparticlc 
bandwidths remains nearly the same and of the order of 
~1.3. However near the Fermi level self-energy effects 
become more pronounced and complicated. In partic- 
ular, the renormalization (bandwidth compression) fac- 
tor grows from 1.3 for stoichiometric composition upto 
nearly 5 for Ko.76Fei.72Se2 (for xz, yz and iz^—r^ bands). 



Though all calculations were done with the same value 
of Coulomb (Hubbard) interaction, this increase of renor- 
malization upon doping tells us that correlations also in- 
crease with doping. Also following recent tendency of ex- 
perimental ARPES papers where energy shifts and and 
renormalization factors are determined separately for dif- 
ferent bare LDA bands [HH to fit ARPES data, we 
have provided quasiparticle mass renormalizations and 
energy shifts, which rather accurately transform bare Fe- 
3d LDA' bands of various symmetries into LDA'-I-DMFT 
quasiparticle bands in different regions of the first Bril- 
louin zone. In fact, the detailed picture of band renor- 
malizations close to the Fermi level in Ki_a;Fe2-ySe2 is 
pretty complicated. For example, the xy,xz,yz bands 
near the F-point and the xz,yz{2) bands in the mid- 
dle of F-X and P-N directions almost do not feel dop- 
ing changes. On the opposite, the xz,yz(l) bands in all 
directions and the 3z^ — r^ band near the F-point be- 
come monotonically more correlated (rcnormalizcd) with 
doping increase. Finally, correlation renormalization of 
the xy band demonstrates non monotonic behaviour with 
doping, becoming more correlated for intermediate dop- 
ings. However, the general conclusion of Ref. [23 re- 
mains valid Ki_3:Fe2-ySe2 systems demonstrates more 
pronounced correlation effects in contrast to 122 iron 



pnictides. This is clearly reflected in the absence of well 
defined quasiparticle bands in the vicinity of the Fermi 
level, which demonstrates a kind of "pseudogap" behav- 
ior ("dark" region around the Fermi level in Fig. 2). The 
interesting problem for the future studies is possible man- 
ifestation of these effects in optical conductivity, as well 
as their role in the formation of superconducting state. 
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FIG. 1: Comparison of LDA'+DMFT calculated spectral 
function maps with renormalized by compression factor of 
1.3 LDA' bands (grey lines) for Ki_iFe2-i;Se2 for different 
hole doping levels: ne=29, 28, 27.2, 26.52 (from bottom to 
top) along high symmetry direction of the first Brillouin zone. 
Fermi level is zero. 



TABLE I: Quasiparticle energy scale renormalization factors and corresponding energy shifts (in eV, in round brackets) for 
different bare Fe-3d LDA' orbitals for all hole doping levels ne in the LDA' scale energy interval from -1.0 eV to 0.4 eV. 



Orbital chracter 


n<i=26.52 


n^=27.20 


ne=28.00 


ne=29.00 


xy 


1.5 (-0.23) 


3.9 (-0.73) 


2.65 (-0.61) 


1.7 (-0.35) 


xz,yz (1) 


4.2 (-0.78) 


3.0 (-0.75) 


2.6 (-0.69) 


1.7 (-0.38) 


xz,yz (2) 


2.3 (-0.48) 


2.5 (-0.60) 


2.6 (-0.69) 


1.7 (-0.38) 


xy,xz,yz 


1.2 (-0.10) 


1.3 (-0.10) 


1.3 (-0.10) 


1.4 (-0.17) 


3z2 - r^ 


4.7 (-0.85) 


2.0 (-0.30) 


1.3 (-0.03) 


1.25 (0.0) 
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FIG. 2: LDA'+DMFT spectral density maps along high sym- 
metry directions of the first Brillouin zone. Maxima of inten- 
sity are shown by stars. Dominant orbital character of quasi- 
particle bands is marked as follows: black squares - xz, yz, 
black circles - xy, white circles - 3z — r , white squares - 
x^ — y^. Fermi level is zero. 
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FIG. 3: Comparison of LDA' (dashed black line) and 
LDA'+DMFT density of states (DOS) for Ki_^Fe2-j;Se2 
for different dopings ne. Thick black line corresponds to 
ne=26.52, thick grey line - ne=27.20, thin black line - 
ne=28.00, grey dot-dashed line - ne=29.00. Fermi level is 




FIG. 4: Bare LDA' (full line) and LDA (dashed line) band 
dispersions for stoichiometric KFe2Se2 system with orbital 
characters explicitly shown. Numbers in brackets after or- 
bital symmetry symbol enumerate corresponding parts of the 
bands in the first Brillouin zone (see Table 1). 



